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Abstract 

Tested morphometries of spiny mice, Acomys, of 181 specimens from 7 localities in Israel and Sinai, to 
evaluate the factors effecting body size differentiation. These specimens represent 3 karyotypes (2 
formal species) comprising 5 populations of Acomys cahirinus (2n = 38); 2 populations of A. cahirinus 
(2n = 36); and 3 populations of Acomys russatus (2n = 66). Each population of the latter is sympatric 
with A. cahirinus. The 7 localities represent a general southward transect of increasing aridity. The 
results indicate that: 1. body characters varied significantly between localities. Body weight and length 
decreased, whereas tail, ear and relative forefoot and hindfoot lengths generally increased with aridity 
in accordance with the Bergmann rule. 2. Morphology was found to be partly correlated with 
allozymic variation. 

The geographic patterns and climatic correlates and predictors of morphological variation in 
Acomys in Israel and Sinai indicate that natural selection, mediated by climatic and biotic factors, is a 
major cause of body size differentiation. The latter in turn contributes to improving thermoregulatory 
efficiency and optimizes the energetics balance. 


Introduction 

Body size is subject to differentiation by both direct and indirect multiple evolutionary 
factors. These involve, among others, ecological (e.g., climatic, habitat, resource availabil¬ 
ity and biotic interactions), physiological (e.g., metabolism, food resources and energe¬ 
tics), demographic (e.g., population size, intra- and interspecific competition) and 
pathogenicity (e.g., parasites and diseases) factors. Their relative importance and complex 
interaction vary spatiotemporally, hence must be evaluated specifically in each case. For a 
general discussion of the functional significance of variation in body size among mammals, 
see Clutton-Brock and Harvey (1983). For a discussion of the general variation and 
correlates of body size in mammals, see Eisenberg (1981). Since heritability of body size 
in mammals is usually high (Shields 1962), geographic variation in size differentiation 
must involve a relatively high genetic component subject to natural selection. The objective 
of this study was to investigate the body size differentiation in spiny mice, genus Acomys, 
in Israel, and hypothesize about its possible determinants. 

Spiny mice, genus Acomys, are tropical murid rodents (Rodentia, Myomorpha, 
Murinae) involving about 18 species ranging in Africa and southwest Asia in rocky 
habitats. Two species of spiny mice, the common spiny mouse, Acomys cahirinus, and the 
golden spiny mouse, Acomys russatus, occur in Israel and Sinai (Shkolnik and Borut 
1969; Haim and Borut 1974, 1975). A. cahirinus is widely distributed in Israel and Sinai 
(Fig. 1), ranging in both mesic and xeric environments comprising the Mediterranean, 
steppe and desert climatic regimes and is thus a climatically euryoek species in range (Table 
1). However, it lives only in rocky areas, and is therefore relatively stenotopic, or narrow 
in niche structure. Two chromosome forms of A. cahirinus, which differ by a single 
Robertsonian change, occur in Israel and Sinai (Wahrman and Goitein 1972, and Fig. 1). 
The Israeli populations possess 2n = 38 chromosomes, and those from Sinai have 2n = 36 
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chromosomes. The two chromosome forms are completely homozygous except for a 
hybrid zone, about 16 km long and 15 km wide, where 2n = 37 hybrids were also found. 
The two chromosome forms are morphologically indistinguishable, but have not as yet 
reached complete reproductive isolation. Although designated as the northern or “Israeli 
Form” and southern or “Sinaitic Form” (Wahrman and Goitein 1972) they may be 
viewed, owing to their chromosomal homozygosity across vast ranges, as derivatives of a 
relatively recent event of speciation, displaying currently its final stages (Nevo 1985). The 
fossil record indicates that A. cahirinus is an upper Palaeolithic colonizer in humid 
Mediterranean Israel, i.e., it presumably appeared some 20.000-30.000 years ago (Tcher- 
Nov 1968). Therefore, its morphological differentiation in Israel is relatively recent and 
traceable over time. 

A. russatus is similar in morphology and life history to the A. cahirinus complex, but it 
differs from the latter in karyotype (2n = 66, Wahrman and Zahavi 1953); in its restricted 
distribution in extreme desert habitats (either hot or cold), in its narrower niche, and in its 
unique physiological adaptations (Shkolnik and Borut 1969; Haim and Borut 1974, 
1975, 1981). It occurs sympatrically with A. cahirinus (either 2n = 38 or 2n = 36) and shares 
the same rocky habitat. However, while A. cahirinus is nocturnal, A. russatus is frequently 
diurnal (Shkolnik 1971). The genetic differentiation and speciation of spiny mice in Israel 
have been recently studied (Nevo 1985) and will be related to the morphometric differenti¬ 
ation in the discussion. Here I present evidence suggesting that body size differentiation in 
Acomys relates primarily to climatic and biotic determinants. 


Material and methods 

Sampling 

A total of 181 (113 males and 68 females) specimens from 7 localities representing 10 populations and 
3 karyotypes were sampled in Israel and in Sinai. These involved 5 (1-5) populations of A. cahirinus 
(2n = 38); 2 populations (6, 7) of A. cahirinus (2n = 36), and 3 populations (8-10) of A. russatus, each 
population being sympatric with a counterpart of A. cahirinus (8 with 4, 9 with 6 and 10 with 7). Data 
on localities and ecogeographical parameters are given in Table 1; distribution and sampling localities 
are shown in Fig. 1. Sampling was conducted in the autumns of 1975 (in Israel) and 1976 (in Sinai) and 
each sample was collected in an area of about one km^. The seven localities of the samples are 
distributed largely along a transect of increasing aridity (see Table 1 and Fig. 1). 


Morphological measurements 

To demonstrate the range of phenotypic variation in the 3 karyotypes of Acomys, 6 measures were 
taken of each adult mouse (and averaged separately for each sex): body weight, and lengths of body, 
tail, forefoot, hindfoot and ear. The ratios of the last four measures to body length were calculated. 


Statistical analysis 

Stepwise multiple regression (SPSS-x 1986) was used to determine whether phenotypic variation of 
morphological traits is influenced by and associated with environmental factors. In addition, 
Pearsonian correlations were computed between all variables including allozymic (see Nevo 1985) and 
morphological variation. The variance in morphology within and between populations was tested by 
analysis of variance, ANOVA. 


Results 

Patterns of morphometric variation 

The means of morphometries of body characters for the 10 populations examined are given 
for males and females separately (Tables 2a and 2b), and the regressions of body weight and 
relative tail and ear sizes on latitude, temperature and the aridity index as expressed by 
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Fig. 1. The distribution of the three karyotypes of Acomys, and the sampled populations: A. cahirinus 
(2n = 38): 1 = Hurfeish; 2 = Beit Oren; 3 = Jerusalem; 4 = Mitzpe Ramon; 5 = Timna. A. cahirinus 
(2n = 36): 6 = Neviot; 7 = Santa Katharina. A. russatus: 8 = Mitzpe Ramon; 9 = Neviot; 10 = Santa 
Katharina 

rainfall are given in Figure 2 in several scatter diagrams. A statistically significant amount 
of variation between localities (in most cases of A. cahirinus, p < 0.001) was found for all 
body characters and ratios (Tables 2a and 2b). The following trends were found in A. 
cahirinus (note that males are better sampled than females). In general, body weight and 
length decreased southward (Fig. 2), whereas the lengths of tail, feet, ear and the ratios of 
all extremities of body length largely increased southward, along with an increase in aridity 
and temperature (Fig. 2). Body weight in A. cahirinus decreases from Hurfeish to Santa 
Katharina by 17.5 % and 19.7 % for males and females, respectively. 

In A. russatus, which is represented by only 3 localities, the following trends were 
found. Body weight and length decreased in both sexes southward from the cool Negev 
highlands of Mitzpe Ramon to the hot Sinai lowlands of Neviot. Body weight decreased 
32.7 % and 27.6 % for males and females, respectively. By contrast, a reverse trend is 























































Table 1. Geographical and climatological data for 7 localities in Israel and Sinai in which 3 karyotypes, Acomys cahirinus^ and Acomys russatus and 

10 populations of Acomys were sampled^ 
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apparent in body weight 
which increases again to¬ 
ward the cold mountain 
tops of southern Sinai in 
Santa Katherina. Body 
weight increased 33.1 % 
and 27.8 % in males and 
females, respectively. 

An opposite trend was 
displayed by the relative 
tail and foot sizes of this 
species. In the better sam¬ 
pled males, they increased 
from Mitzpe Ramon to 
hot Neviot, but decreased 
again in cold Santa 
Katharina (see trends in 
Table 2). These reverse 
trends from Neviot to 
Santa Katharina are incon¬ 
sistent in A. cahirinus. 


Environmental correlates 
and predictors of mor¬ 
phological variation 

Pearsonian correlations 


In the search for environ- | 
mental correlates of mor¬ 
phological variation we 
used several geographic 
and climatic variables 
(Table 1): Longitude = 
Ln; Latitude = Lt; and Al¬ 
titude = Al, were used as 
geographic variables. 
Temperature (mean annual 
= Tm; mean January = Tj; 
and mean August = Ta) 
and moisture (mean annu¬ 
al rainfall = Rn; mean mid¬ 
day relative humidity = 
Hu; and mean annual 
evaporation = Ev) were 
used as climatic variables. 
To assess the ecological 
background structure the 
correlation matrix of the 
ecological variables was 
analysed. Following are 
some correlations, where 






Table 2a. Morphometries of adult Acomys cahirinus and Acomys russatus accompanied by analysis of variance (ANOVA); Males 
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Table 2b. Morphometries of Acomys cahirinus and Acomys russanus accompanied by analysis of variance (ANOVA); Females 
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LATITUDE AUGUST TEMPERATURE ANNUAL RAINFALL 

Fig. 2, Regressions of body weight and of relative ear and tail size of Acomys cahirinus and A. russatus 
on latitude, August temperature and annual rainfall, a Population mean of males ai Acomys cahirinusy 
o population mean of females of Acomys cahirinuSy a population mean of males of Acomys russatusy 
• population mean of females of Acomys russatuSy □ both sexes of Acomys cahirinuSy (M) males, 
(F) females 


figures in parentheses are Pearsonian r’s for all 7 localities of Acomys. All 3 temperature 
variables, Tm, Tj and Ta, are highly correlated (0.88-0.97); A1 and Tm (-0.97); Rn and Ta 
(-0.29); Rn and Hu (0.92); Ev and Ta, Rn and Hu (0.55, -0.86, -0.94, respectively). The 
low correlation between Rn and temperature is due to the reversed temperature trend in 
Santa Katharina. When the latter locality is excluded, the correlation between Rn and 
temperature is negatively high as usual across the studied area (r between Rn and the 
temperature variables are from -0.67 to -0.86). 

In general, rainfall increases and temperature decreases northward with latitude; r 
(Rn-Lt) = 0.94; N = 7; and r (Tj-Lt) = -0.93 for 6 localities excluding Santa Katharina. In 
other words, aridity increases southward toward the Negev and Sinai deserts. Tm and Hu 
were eliminated from the multiple regression analysis in order to reduce the level of 
intercorrelations between the environmental variables. 
































Table 3. Pearson correlation between morphology and climatic variables in Acomys cahirinus 7 populations 
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Morphological correlates with climate 

As can clearly be seen in Table 3, high and 
significant correlations, above those expected 
by chance, relate morphometries and climatic 
factors. Body size significantly decreases (r, 
for males, = -0.79; p < 0.05) southward with 
evaporation whereas the relative tail size in 
males increases southward with evaporation 
(r = 0.96; p < 0.001). Even single climatic 
factors explain significantly the morphologi¬ 
cal variance in body ratios, but much more of 
the variance is explained by two or three 
variable combinations (see later, multiple re¬ 
gression analysis. Table 5). 

Morphological correlates with allozymes 

The frequencies of some allozymes show sig¬ 
nificant correlations with morphometric 
means (Table 4). The analysis is based on the 
assumption that there is no real difference in 
allozyme frequencies between the sexes. The 
reliance on male morphometries stems from 
their better sampling. Because of the high 
correlation between the morphometries of 
both sexes, the analysis of females gave simi¬ 
lar results. The following alleles showed sig¬ 
nificant correlations with morphometries, 
the numbers in parentheses indicate the 
number of body variables correlated with the 
alleles: Aat-1^ (3); Est-2^ (2); Est-2^ (1); Est- 
3" (1); Est-5^ (3); Est-b"" (1) and Est-Z’^ (2). In 
addition, in females only: Pept-2‘^ (r = -O.SE*' 
= p<0.05 with forefoot); Est-b*^ (r = 0.78''*' 
with ear, and r = -0.82'*' with hindfoot-body 
ratio). It thus appears that allozymic and 
morphological diversities are at least partly 
significantly correlated, primarily with the 
allele frequency in several esterase loci. Both 
allozymic and morphological variations are 
correlated with environment. 

Multiple regression analysis 

A test of the best predictors of the mor¬ 
phological variables of the populations of the 
A. cahirinus complex was conducted by step¬ 
wise multiple regression analysis, MR (SPSS- 
X 1986), employing the aforementioned var¬ 
iables as dependent variables and ecogeo- 
graphic factors as independent variables. The 
results are given in Table 5. A substantial 
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Table 5. Coefficients of multiple determination (R^) employing as dependent variables 6 morpho¬ 
logical variables and 4 morphological ratios of Acomys cahirinus and as dependent variables 
geographic + climatic, and climatic factors 


Variable 

Stepwise model 

Geographic + Climatic variables Climatic variables 

Weight 

Same as for climatic variables 

Ev 

EvRn 

EvRnTa 

dd 


0.53^ 

0.73^ 

0.84 

Body length 

Ev EvLn EvLnAl 

Ev 

EvRn 

EvRnTj 

dd 

0.62'*- 0.72^ 0.95"' 

0.62’=' 

0.68^ 

0.70 

Tail length 

Same as for climatic variables 

Ev 

EvRn 

EvRnTa 

dd 


0.85’='’=' 

0.90’=' 

0.93’=' 

Forefoot length 

Ta TaLn TaLnLt 

Ta 

TaEv 

TaEvRn 

dd 

0.39 0.54 0.85^ 

0.39 

0.43 

0.45 

Hindfoot length 

Ln LnRn LnRnEv 

Tj 

TjTa 

TjTaRn 

dd 

0.24 0.38 0.83 

0.20 

0.64 

0.65 

Ear length 

Same as for climatic variables 

Rn 

RnTa 

RnTaEv 

dd 


0.80’-"- 

0.90’=' 

0.91’=' 

Tail/Body ratio 

Ev EvLn EvLnAl 

Ev 

EvTa 

EvTaTj 

dd 

0.92''^'’=- 0.96-’^’^ 0.98^'^' 

0.92’•'’='’=' 

0.93’=-’=' 

0.97’=' 

Forefoot/Body ratio 

Ev EvLt EvLtAl 

EV 

EvRn 

EvRnTa 

dd 

0.75"' 0.93"''-' 0.96"' 

0.75’=' 

0.87’=' 

0.91’=' 

Hindfoot/Body ratio 

Same as for climatic variables 

Ev 

EvRn 

EvRnTj 

0.94’=' 

dd 


0.58’=' 

0.80’=' 

Ear/Body ratio 

Rn RnAl RnAlTa 

Rn 

RnTj 

RnTjTa 

dd 

0.87’='’=' 0.90’='-'^ 0.94’=' 

0.87’=""' 

0.90’=-’=' 

0.91"' 

Abbrev. of environmental variables: Ln = Longitude, Lt = 

= Latitude, A1 

= Altitude, Ta = mean 

August temperature; Tj 
evaporation. ^ = p<0.10, 

= mean January temperature, 

’=' = p <0.05, '"""' = p < 0.01, ’='’•'’=' 

Rn = annual 
= p <0.001. 

rainfall, 

Ev = annual 


amount of the geographic variation in body characters was explained by the climatic factors 
of water availability (evaporation and rainfall) and temperature. Note that the variance in 
body ratios, which represents the relative estimates of body extremities (tail, ear, fore- and 
hindfoot lengths) is significantly explained by a single factor (either evaporation or 
rainfall), but a two or three variable combination explains above 90 % of the variance. 


Discussion 

I will subdivide the discussion into tow parts, 1. general, including other examples of 
geographic variation of body size in both space and time across Israel and in other regions, 
and 2. specific, analysing the results of the Acomys complex. 

General comments on body size and energetics 

The usual explanation of the positive correlation of body weight with latitude in 
homoiotherms (Bergmann’s rule) relates to the physics of heat exchange, providing a 
thermoregulatory device for conserving energy in colder climates by a larger body size, 
and dissipating energy in warmer climates by a smaller size. This interpretation of 
Bergmann’s rule was criticized by Scholander (1955) on several grounds. First, many 
species did not follow the rule and those that did show clinal increases in weight were 
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physiologically unimportant. Second, many of the characters supposedly associated with 
the rule have no significance in heat exchange. The rule was also questioned when some 
poikilotherms reflected it (Ray 1960). Rensch (1932) concluded that when factors other 
than temperature affect body size the maximum is reached in the optimal portion of the 
species range. Mayr (1970) emphasized that diametrically opposed size trends occur in 
species differing in life cycle patterns, sexual maturity, and optimal habitats. 

In a series of studies on the patterns and evolution of endothermy in the phylogeny of 
mammals and other vertebrates McNab (1970, 1971, 1974, 1979, 1980, 1983) has critically 
analysed the influence of body size on the energetics, food habits and population biology 
of mammals. In his analysis of the ecological significance of Bergmann’s rule, McNab 
(1971) emphasized that those species which do conform to the Bergmann rule “are usually 
carnivores or granivores; a change in their body size reflects a change in the size of their 
prey. A latitudinal change in the size of available prey is due either to the distribution of 
the prey species or to the distribution of other predators utilizing the same prey species. 
Only the smallest species of a set of similar predators normally will conform to Bergmann’s 
rule, and then only beyond the limits of distribution of the largest species. These changes 
in size seem to be another example of character displacement.” 

McNab concluded that Bergmann’s rule usually is a special case of a more general 
phenomenon related to food frequency and its size and the presence of other species that 
utilize the same food resources. In a recent analysis McNab (1983) analysed the complex 
relationship between energetics and body size and the limits to endothermy. Most 
vertebrates follow a Kleiber relation down to a “critical” mass, below which the scaling of 
metabolism must be changed to ensure the maintenance of endothermy. Critical mass 
varies inversely with the level of energy expenditure. The perfection of endothermy may 
always require an evolutionary decrease in mass. These considerations may not only be 
true at higher taxonomic levels, but also intraspecifically. 

Body size as a “niche difference” 

Body size variation may permit the coexistence of close species, thus providing an 
important niche parameter. The reason is that different sized animals eat different sized 
foods, hence utilize differential resources, at least partly nonoverlapping in size, to allow 
coexistence (MacArthur 1972). Wilson (1975) describes patterns of convergence and 
divergence of body size unifying in his model seemingly contradictory concepts: 1. the 
concept that differences in body size promote a niche difference (Hutchinson 1959; 
Brown and Wilson 1956; Schoener 1965, 1967, 1970, 1974a, b; Grant 1968, 1972), and 
2. the concept that differences in body size set up a competitive gradient whereby the larger 
can exclude the smaller (Brooks and Dodson 1965; Galbraith 1967). 

Spatiotemporal variation in size of other rodents in Israel 

Several studies analyse the spatial and temporal morphological variation in Israeli rodents. 
Size variation of Meriones tristrami was analysed by Chetboun and Tchernov (1983). Its 
size variation over the geographic range of Israel and northern Sinai and during the late 
Middle and Upper Quaternary of Israel is correlated with environmental factors. The 
allopatric fossil populations were found to fluctuate in accordance with Bergmann’s rule 
(i.e., the maximum body size was obtained during colder phases). Recent populations did 
not show north-south body size gradient of high correlation with climatic factors. A 
sudden shift in body size of M. tristrami is observed when sympatry with its congener M, 
sacramenti occurs in the Coastal Plain, and another (and smaller) shift is noticeable when 
sympatry with a third congener (M. crassus) takes place in the Negev region. The allopatric 
convergence and sympatric divergence suggest an ecological interaction (character dis- 
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placement) among the congener species. A significant north-south increase in the relative 
volume of the bulla tympanica of M. tristrami is shown, with the mastoid portion playing 
the main role in geographic variation. The bulla volume is highly correlated with climatic 
factors and shows no divergent shifts within the sympatric region. It is suggested that in 
fossil populations the bulla tympanica can be used for evaluating palaeoclimatic fluctua¬ 
tions. 

Size trends and Pleistocene paleoclimatic response of the murid genus Apodemus 
(subgenus Sylvaemus) were studied by Tchernov (1979). Eleven Pleistocene sites in Israel 
yielded samples of successive populations of Apodemus. Three species occur: two {A. mys- 
tacinus and A. sylvaticus) are continuously represented in the area, while the third, 
A. flavicollisy occurs only during colder periods. The separation of fossil populations is 
possible only on biometrical grounds. The observed distribution of coexisting populations 
revealed significantly different modes of each parameter, caused by strict ecological 
exclusion. Molars showed three well-defined morphological size-groups within each of the 
three species when they coexisted in closely related adaptive zones. Neither significant 
complication, simplification nor clinal variation occurred in tooth structure since early 
Middle Pleistocene. 

In the Middle East three different lineages of Apodemus, which never overlapped in 
size, were present during the Middle and Upper Pleistocene. When all three species 
coexisted, parallel changes in size with time occurred in all of them simultaneously, 
apparently in accordance with Bergmann’s rule. Yet when Apodemus flavicollis was 
temporarily eliminated from the scenario in warm periods, the other two species filled the 
vacated space independently of climate, and in contrast to Bergmann’s rule. Two factors 
greatly affected the character displacement of the three chrono-species: interspecific 
competition and climate. As competitive exclusion of these three closely related species is 
and was heavy, it is suggested that no factor other than size was involved in their evolution 
during the last 700,000 years. The significant size changes of each of these lineages, with no 
overlap among them, were interpreted by correlations with climatic factors and by 
ecological exclusions. 

In an extensive survey of faunal turnover and extinction rate in the Levant, Tchernov 
(1984) analysed size changes in the late Pleistocene and postglacial period. Size change is 
exemplified most drastically by postglacial radical dwarfing seen in several species of 
mammals in Israel: Canis lupus, Ursus arctos, Sus scrofa, and Felis sylvestris (Kurten 1965); 
Spalax ehrenhergi and Microtus guentheri (Tchernov 1968); Gazella gazella, Vulpes 
vulpes, and Canis lupus (Davis 1977). Davis showed that the Aurignacian wolf of Israel 
attained the same dimensions as the present population of southern Sweden. Differences in 
mean temperature (for the hottest or coldest months) between the two regions is around 
15 °C. According to Davis (1977), dwarfing commenced within the Natufian period 
(10,000 to 12,000 B.P.) and was completed approximately during the Neolithic. It is 
obvious that a conspicuous climatic change took place in the eastern Mediterranean region, 
an effect not clearly shown by the mammalian fauna turnover. 

A detailed analysis on differentiation of body size in general (Nevo et al. 1986) and on 
42 skull and body characters (Nevo et al. 1988) have been conducted on the four 
chromosomal species (2n = 52, 54, 58 and 60) of subterranean mole rats of the Spalax 
ehrenhergi complex in Israel. For differentiation of body size, the weight of 1,653 
subterranean mole rats comprising 12 populations of the four chromosomal species were 
analysed. The results indicated that there is a southward latitudinal gradient in body size. 
Northern animals living in cooler and more productive mesic environments are larger than 
southern animals living in warmer and less productive xeric environments. Body size is 
negatively correlated with temperature variables, and positively correlated with plant cover 
(reflecting productivity or food resources) and rainy days. The best predictors of body 
size, explaining up to 87 % of the variation in size included various combinations of 


Natural selection of body size differentiation in Acomys 


93 


temperature variables and plant cover. The conclusions were that in both adaptation and 
speciation natural selection is a major agent of differentiation of body size in accordance 
with multiple factors, primarily temperature and food resources operating on the energe¬ 
tics balance. 

In a complementary study (Nevo et al. 1988) the morphometries of 327 adult 
subterranean mole rats comprising 44 populations distributed across the ranges of the four 
chromosomal species of the Spalax ehrenhergi superspecies in Israel has been studied on a 
total of 42 skull and body variables. The results indicated that climatic variable combina¬ 
tions of temperature and water availability explain a significant part of the variance of most 
skull and body variables. Morphological diversity displays a southward gradient of 
decreasing size in skull and body variables which are significantly explained by climatic 
selection (and indirectly by decreasing resource availability southward). Hence it is 
adaptive and explicable on even very low selective pressures over evolutionary time. 

Geographic variation of decreasing body size southward occurs in Israel in other 
homoiotherm mammals (e.g., wolves, leopards and hedgehogs; in Mendelssohn 1982) 
and in birds (partridges; in Nisani 1974). In contrast, increase in body size southward 
occurs in some poikilotherms (e.g., toads; in Nevo 1972). The decrease in body size 
southward in mammals and birds may be due to the sole or combined effects of higher 
temperatures and lower productivity in the xeric southern Negev and Sinai deserts. The 
increase in size of toads southward may be due to the fact that larger toads, having a 
relatively smaller evaporative body surface area, are capable of withstanding longer periods 
of dessication and are therefore selectively superior in arid habitats. 


Ecophysiological, genetical and morphological differentiation of spiny mice, Acomys 

What are the evolutionary forces molding body size differentiation in spiny mice in Israel? 
I will first summarize the ecophysiological background structure, then discuss body size 
differentiation, relate it to allozymic diversity and assess the evidence indicating that 
natural selection is a major architect in population differentiation at both the molecular 
(i.e., genetic) and organismal (i.e., morphological) levels. See Endler (1986) for a critical 
review of natural selection in the wild. 


Ecophysiological background 

A. cahirinus is a climatically generalist species ranging over widespread mesic and xeric 
habitats across Israel and part of the Sinai, whereas A. russatus is a habitat specialist species 
restricted to extreme xeric desert conditions in these regions. Recently, A. russatus was 
described from the colder mountain tops of the southern Sinai desert at altitudes of up to 
2650 m (Haim and Tchernov 1974). Physiological analysis of the Sinai mountain top 
populations of A. russatus indicates that it is adapted to cold-stress physiology by means of 
a higher nonshivering thermogenesis (Haim and Borut 1974, 1975, 1981), whereas 
animals from the En-Gedi Judean Desert population were unable to produce enough heat 
to keep warm below 18 °C (Shkolnik and Borut 1969), However, both species {A. 
cahirinus and A. russatus) are narrow-niched species restricted to life in rocky crevices 
where microclimatic conditions in their microhabitats are relatively constant and far more 
moderate than the surrounding climates (Schmidt-Nielsen 1964; Shkolnik 1966); yet 
they differ in their thermal niche-breadth (see below). 

2. A. cahirinus is able to regulate its body temperature over a much greater range of 
ambient temperatures than A. russatus in accordance with its larger habitat range (Weis- 
SENBERG 1977; Weissenberg and Shkolnik 1977). However, A. russatus is physiologi¬ 
cally better adapted for a diurnal life than nocturnal A. cahirinus by its ability to withstand 
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high ambient temperatures and generate less heat by metabolism, as well as its capacity to 
derive water from salty, succulent plants (Shkolnik and Borut 1969; Shkolnik 1971). 

3. Both species are adapted to desert conditions through an adaptive syndrome which is 
more pronounced mA. russatus than in A. cahirinus. Their metabolic rates deviate from the 
expected based on weight, 34.5 % in the former and 13.5 % in the latter, thereby reducing 
the water required for evaporation. Likewise, maximum chloride concentration is higher in 
the urine of A. russatus (1500 mN), and it is almost twice that of A. cahirinus. Finally, both 
species have among the very highest mammalian urine concentration (4700-4800 mN; 
Shkolnik and Borut 1969). 

4. High ratios of evaporative water loss characterize both species, but more so in 
russatus than in A. cahirinus., involving a high cutaneous water loss, about 60-70 % of the 
total. The high cutaneous water evaporation, particularly in diurnal A. russatus, allows it to 
dissipate at an ambient temperature of 30 °C, over the third of the heat it generates. This 
efficient cooling device through a high evaporation rate, which appears very high even 
when compared to other diurnal rodents, is compensated by a very efficient kidney, by 
food choice comprising snails and succulent plants and by behavioral adaptations of shade- 
path selection (Shkolnik 1971). All four ecophysiological factors suggest a tropical (hot 
and humid) evolutionary origin of Acomys (Shkolnik 1966), as is also substantiated by 
taxonomic (Ellerman and Morrison-Scott (1951), paleontological (Tchernov 1968, 
1975, 1984) and cytogenetic (Warhman and Goitein 1972) studies. 


Genetic and morphological patterns 

The genetic differentiation and speciation in Acomys has been described elsewhere (Nevo 
1985). Here I summarize the major genetic and morphological patterns. 1. Low levels of 
genic diversity (A, P, H) characterized Acomys but they were higher in the A. cahirinus 
complex than in A. russatus. 2. A substantial proportion of the variant loci, 40-50 % were 
either localized or sporadic, suggesting sharp local and regional differentiation of alleles 
across the mesic-xeric ranges of Acomys in Israel and Sinai. 3. In both A cahirinus and A. 
russatus the level of polymorphism and number of alleles increased southward with aridity. 
4. Levels of A, P and H varied at enzyme and protein loci in different functional classes. 
Loci whose enzymes utilize substrates originating from the external environment are far 
more genetically diverse than loci whose enzymes utilize internal metabolites (Gillespie 
and Kojima 1968; Johnson 1974). Likewise, regulatory enzymes were more variable than 
nonregulatory enzymes Qohnson 1974). 5. Significant linkage disequilibria were found in 
A. cahirinus in three populations involving six alleles. 6. Deviations from Hardy-Weinberg 
equilibria owing to heterozygote paucity were found in eight populations involving five 
loci. 7. Significant heterogeneity between 14 polymorphic loci in their effective inbreeding 
coefficients suggest the operation of natural selection. This test was criticised on statistical 
grounds (Ewens and Feldman 1976) and I use it here only as a supportive, and not a 
conclusive result. 8. Polymorphism and allozymic variant at six polymorphic loci were 
correlated with and predicted by climatic factors of water availability and temperature. 9. 
Body characters vary significantly between geographical localities: Body weight and length 
decreased, whereas tail, ear and relative fore- and hindfoot lengths increased and were 
correlated significantly with aridity and temperature. 10. Morphology was found to be 
partly correlated with allozymic variation. 11. Genetic patterns of A. cahirinus and A. 
russatus varied in three sympatric localities. 12. Mean genetic distance, D, was very small 
within A. cahirinus (2n = 38 and 2n = 36), and within A. russatus. It was also very small 
between the karyotypes of A. cahirinus (2n = 38 and 2n = 36), but high between the A. 
cahirinus complex and A. russatus. 
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Adaptive differentiation of body size in spiny mice 

The interspecific size differences between nocturnal A. cahirinus and diurnal A. russatus 
may relate to their striking pattern of competitive exclusion (Shkolnik 1966, 1971) which 
results in opposed patterns of circadian activity. This variation of body size between the 
species may represent niche differences (see Wilson 1975). The heavier body weight of A. 
russatus^ as compared with A. cahirinus (Table 2) may be associated with its diurnal life, 
where a relatively small surface to volume ratio may be adaptively superior in dry and 
high-radiation environments. 

In general, in both A. cahirinus and A. russatus^ sizes decrease and extremities increase 
with aridity and temperature (Fig. 2). These two trends reflect a positive correlation of 
weight with latitude (Bergmann’s rule) and a negative correlation of body extremities with 
latitude (Allen’s rule). Both patterns appear to better adapt Acomys to a progressively 
increasing heat load southward by improving its thermoregulation capacity through larger 
dissipating heat surfaces. Most fossorial mammals conform to these trends, although the 
interpretation of the trend in other mammals is often complex, and may involve inter¬ 
species interaction or character displacement (e.g. McNab 1979; Chetbourn and Tcher- 
Nov 1983; Tchernov 1979; Wilson 1975). 

Although one could suggest character displacement as a factor involved in body size 
differentiation of spiny mice, it appears very unlikely. In all three populations where both 
A. cahirinus and A. russatus coexist (localities 4-8, Mitzpe Ramon; 5-9, Neviot; and 7-10, 
Santa Katharina) the size difference appears to reflect an extension of general regional 
gradient of each of the species, rather than local character displacement, as is evident by the 
general ranking pattern in each of the species and in both sexes, i.e., decrease in size 
southward from the Central Negev desert (Mitzpe Ramon) to the Sinai desert (Neviot). 
The increase in size of A. russatus in both sexes at the mountaintop of Sinai (Santa 
Katharina) appears to support the thermoregulatory hypothesis. The contribution of a 
bigger body size to thermoregulatory efficiency in mountaintop colder environments 
complements the physiological adaptation by nonshivering thermogenesis of this popula¬ 
tion (Haim and Borut 1974, 1975, 1981). 

The involvement of natural selection in genetic and morphological population differ¬ 
entiation is suggested by points number 1-11 mentioned above for both allozymes and 
body characters, which are intercorrelated and are explicable at least partly by the 
environment. Neither gene flow nor genetic drift can explain satisfactorily the genetic and 
morphological patterns found. If gene flow was a major factor, sporadic and localized 
alleles would not comprise about half of the variant alleles, and dines should have been 
more abundant. If genetic drift was a major factor, fixation of alternative alleles should 
have been more pronounced within the A. cahirinus complex and allele distribution would 
not be explained, at least partly, by climatic factors. Likewise, it is unlikely that small size 
population effects are responsible for the regional patterns. Both species are very abundant 
in the continuous rocky habitats across Israel and parts of Sinai. The intercorrelation 
between allozymic and morphological diversities reinforce their adaptive component. 

In sum the morphological differentiation of body size, like the genetic one (Nevo 1985) 
appears to relate to environmental factors. Climatic factors appear substantial in molding 
body size as a contributor to thermoregulatory efficiency interacting with genetic, 
physiological and behavioral factors to optimize the energetics balance. Additional factors 
such as decline in resource availability southward may also contribute to the decrease in 
body size toward the hot Negev and Sinai deserts. It is noteworthy that desert rodents in 
Israel show a peaked curve of diversity over productivity, drastically declining in the deep 
desert (Abramsky and Rosenzweig 1984). 

While in general temperature and aridity (i.e., lower productivity) increase southward 
across the Acomys range, they are correlated and impossible to disentangle. Partial 
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separation of both factors may be achieved between hot and arid Neviot populations and 
cold and arid Santa Katharina populations, both in the Sinai desert. Here, while aridity 
(hence productivity) is similar, temperature differences are extreme. Therefore, size 
differentiation between both sites appear largely to be determined by temperature, in 
accordance with the Bergmann rule. This is observed primarily in A. russatus. However, 
the relative importance of declining food resources in the southward decrease in size 
remains to be further evaluated by critical, specifically designed experiments. 

Natural selection appears to be a major architect of size as well as allozyme (Nevo 1985) 
differentiation in spiny mice and of their interaction. A critical discussion of the operation 
of natural selection in natural populations has been presented by Manly (1985) and 
Endler (1986). 

Evolutionary history and the rate of morphological differentiation 

The evolutionary center of origin of spiny mice based on taxonomic grounds (Ellerman 
and Morrison-Scott 1951; Wahrman and Goitein 1972 and references therein) is the 
Ethiopian Region, and Israel represents a northern extension from this center. The African 
genus Acomys existed apparently in the Near East since the Pliocene, as it is found in 
Cyprus which separated from the mainland not later than the early pleistocene. This 
Pliocene ancestor presumably radiated into several closely related species (e.g. Atallah 
1967). No Acomys remains were found however, in the Pleistocene of Israel until as late as 
the Aurignacian in the upper Palaeolithic, i.e., some 20,000 years ago (Tchernov 1968, 
1975). A. cahirinus is a recent colonizer to the mesic but deforested Mediterranean region, 
particularly in the Epipaleolithic period, i.e., 15,000 B.P. Noteworthy, in the Mediterra¬ 
nean region A. cahirinus occurs primarily on warm and dry south-facing wadi slopes, 
presumably reflecting its derivation from steppic or desert origins, whereas the cooler and 
more humid north-facing slopes in Mount Carmel, for example, are primarily inhabited by 
the European Apodemus (Nevo, unpubl.). A. russatus may have originated in the extreme 
southwest Asiatic arid zones. 

The three chromosomal forms of Acomys in Israel, A. cahirinus (2n = 38, 36) and 
A. russatus (2n = 66) appear to represent two speciation events: an old and a recent one. 
This deduction is based on both karyotypic differentiation (Wahrman and Zahavi 1953) 
and genetic distances (Nevo 1985). The genetic distance separating A. russatus and A. 
cahirinus is D = 0.30. A rough evolutionary divergence time for a pair of species estimated 
from electrophoretic data based on their genetic distance can be obtained by t = 5 x lO^D 
(Nei 1975). On this basis, A. russatus and A. cahirinus separated 1,500,000 ± 50,000 years 
ago. On the same basis, and in sharp contrast to this older speciation event, the karyotypic 
differentiation within A. cahirinus occurred 115,000 ± 40,000 years ago. This figure is in 
accordance with the fossil record which reveals A. cahirinus in Israel only as late as 20,000 
years ago (Tchernov 1975). It supports the hypothesis that the Israeli form of A. cahirinus 
(2n = 38) has probably been derived from the Sinai tic form (2n = 36) through fission of 
one metacentric chromosome, rather than the reversed, cytologically simpler process 
(Wahrman and Goitein 1972). 

The evolutionary recency of the Acomys cahirinus complex in Israel, i.e., 20,000 years, 
suggests that natural selection must have been instrumental in driving a relatively high 
evolutionary rate of size differentiation. 
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Zusammenfassung 

Naturliche Selektion in der Differenzierung von Kdrpermajlen hei der Stachelmaus, Acomys 

Untersucht und miteinander verglichen wurden 181 Stachelmause des Genus Acomys von 7 Lokalita- 
ten aus Israel und dem Sinai. Diese Tiere entsprechen 3 Karyotypen (2 Arten), welche 5 Populationen 
von Acomys cahirinus (2n = 38), 2 Populationen von A. cahirinus (2n = 36) und 3 Populationen von A. 
russatus (2n = 66) umfassen. Jede Population von A. russatus ist mit einer von A. cahirinus 
sympatrisch. Die 7 Lokalitaten erstrecken sich von Norden nach Siiden und entsprechen Regionen mit 
schnittweise zunehmender Ariditat. Die Ergebnisse der Untersuchungen zeigen: 1. Die Tiere der 
verschiedenen Lokalitaten unterscheiden sich in den erfa£ten korperlichen Merkmalen signifikant 
voneinander. In Ubereinstimmung mit der Bergmann’schen Regel zeigen Korpergewichte und -langen 
Abnahmen mit zunehmender Ariditat, die absoluten Langen von Schwanz und Ohr sowie die 
relativen von Vorder- und Hinterfuft nehmen hingegen allgemein zu. 2. Die morphologischen 
Besonderheiten sind teilweise korreliert mit enzymatischen Variationen. Die Beziehungen zwischen 
geographischen Faktoren sowie klimatischen Besonderheiten und morphologischen Unterschieden 
bei Acomys in Israel und dem Sinai zeigen an, daB die naturliche Selektion als wichtigster Faktor fiir 
die Differenzierungen der K6rperma£e angesehen werden mu£. Diese Differenzierungen tragen bei 
zur Verbesserung der Thermoregulation. Durch sie wird die Energiebalance optimiert. 
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